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Abstract 

A theoretical study of the structural transition of impurity nanorod array in epitaxial 
YBa2Cu307_5 films on vicinal SrTiOa substrates is presented. Two possible types of film/substrate 
interface were considered with one assuming a complete coherence while the other, defective as man- 
ifested in presence of anti-phase grain boundaries. Only in the former case, the increase of the 
vicinal angle of the substrate leads to a substantial change of the strain field in the film, resulting in 
a transition of the nanorod orientation from the normal to in-plane direction of the film. The cal- 
culation of the threshold vicinal angle for the onset of the transition and the lattice deformation of 
the film due to the inclusion of the impurity nanorods is in very good agreement with experimental 
observations. 
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Nanocomposites, in particular those with the nanostructured second phase doped into a 
primary matrix, have received much attention recently since functional performances of the 
material may be directly enhanced through microstructure engineering. Examples of such 
nanocomposites include ferroelectric and multiferroic oxides, magnetic oxides, and supercon- 
ducting cupratesi"-. For high-temperature superconducting YBa2Cu307_5 (YBCO) films, 
many experimental efforts have been made to generate impurity columnar nanostructures 
in the film to maximize the superconducting critical current (Jc) in magnetic fields. The 
doping of insulating impurities such as BaZrOs (BZO), BaSnOs (BSO), YBa2Nb06, and 
rare-earth tantalates have been identified as promising possibilities for that purpose^i^*^"-. 
While extended columnar defects provide an increased overall pinning force in the supercon- 
ducting films as they effectively pin magnetic vortices along their length, a splay alignment 
of the columnar defects has been found to reduce the possibility of vortex hopping and re- 
sult in stronger pinning in high magnetic fields^i"— . The tuning of the defect alignment in 
YBCO films could therefore provide a possibility to optimize the overall behavior of Jc with, 
especially, reduced Jc anisotropy with respect to the direction of magnetic field. In order 
to achieve the controlled alignment of columnar nanostructures such as BZO nanorods in 
YBCO films, we recently employed vicinal SrTiOa (STO) substrates to introduce control- 
lable structural parameters for the formation of nanostructuresi^*^. It was observed that 
the orientation of BZO nanorod arrays in YBCO films depends on the vicinal angle of the 
vicinal substrate. Increasing the vicinal angle can change the nanorod alignment from the 
normal direction of the film, to splayed around the normal at 5°, and to the ah plane of 
the film at larger angles^*^. In the splayed case, much enhanced overall Jc was obtained 
in magnetic fields up to 5 T with almost any field orientation. Theoretical understanding 
of the role of the vicinal substrate on the nanorod orientation is thus important but not 
currently available. This paper presents the first theoretical study on the impurity nanorod 
alignment in epitaxial YBCO films on vicinal substrates. The propose model for the for- 
mation of nanostructures is rather generic and could be applied to a variety of epitaxial 
nanocomposite films to guide growth of nanostructures in a controllable fashion. 

In this study, the effect of vicinal substrates on the formation of impurity nanorods in 
epitaxial YBCO films is modeled in terms of the effect of lattice-mismatched substrates. The 
configurations of the nanorods in epitaxial YBCO films on lattice-mismatched substrates 
have been studied with a micromechanical modeU^, in which the formation of the nanorods 

2 



was assumed to be the consequence of the relaxation to the energetically-preferred elastic 
equilibrium of coherently strained lattices due to lattice mismatches among film, dopant, 
and substrate. The possibility to fabricate impurity nanorod arrays with the vertical or 
horizontal alignment in the c-oriented YBCO films depends on the elastic energy of the 
strained lattices with respect to other possible nanorod configurations. The elastic energy 
of strained lattices is defined as 

Eei = J EidV+ J E2dV+ J EsdV (1) 

film dopant substrate 

where Ei is the elastic energy density for film {i = 1), dopant {i = 2), and substrate (i = 3), 
respectively. Considering the tetragonal symmetry of twinned YBCO film and the cubic 
symmetry of substrate and dopant, the elastic energy density can be written as 

Ei = ^AjiMn + ^Ai2M33 + Ai3MiiM33 + Ai4M?3 (2) 

where Ujk is the strain tensor and = Cil, Xi2 = C33, Ai3 = c^g, and Aj4 = C55 are the 
elastic constants of the material labeled with i. The elastic constants used in this study can 
be found in Refs.— . Note that the interaction energies at film/substrate and film/dopant 
interfaces are included implicitly in Ei as Ujk is the solution of equilibrium equations with 
the boundary conditions that are the result of the interface interactions. The equilibrium 
equations can be written as 

where Xi is the coordinate along the [100] {i = 1), [010] {i = 2), and [001] [i = 3) direction, 
respectively, and ajk = dEi/dujk is the stress tensor. At an interface between two coherently 
bonded lattices, the boundary condition of Eq. ([3]) prescribes continuity of the force on the 
interface and allows for a discontinuity of the strain across the interface, i.e.— 



k 

5]sfe[a,fc(l) + (T,fc(2)] = (4) 
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where n and s are the unit vectors normal and tangential to the interface, {ujk{l),(Jjk{l)) 
and {ujk(2),ajk{2)) are the strain and stress at the interface in lattice 1 and 2, respectively, 
and fk is the lattice mismatches at the interface along the Sk direction. At the top surface of 



a film the boundary condition is simply UkCTjk = and deep inside the substrate the strain 
vanishes. Consider the case that the c-oriented YBCO film mismatches with the substrate 
lattice along the [100] direction. To determine the energetically-preferred orientation of 
impurity nanorods in the film, the equilibrium strain was solved from Eq. for the 
nanorods aligned in the [001] or [100] direction of the film and the difference of the elastic 
energies with the two alignments were obtained as^ 



p+—p)G ^ 



(5) 



where i?ez [100] and -Eg; [001] are the elastic energies of the film with the nanorods aligned in 
the [100] and [001] directions, respectively, Eq is a positive function, p is the volume density 
of nanorods, 

= Ail - Ay Ai2 (6) 

is of the elastic constants of the film {i = 1), dopant {i = 2), and substrate {i = 3), 
respectively, and a = A12/A11. In Eq. ([5]), 
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a, (7) 



{Wi +W2){f3 + Sf3) 

where /i + 6fi and /s + 5/3 are the lattice mismatches between the dopant and the strained 
film lattice due to the mismatched substrate. These lattice mismatches were calculated in 
Ref.^ as 

where (ai,Ci) and (025 C2) are the lattice constants of the film and dopent along [100] and 
[001] direction, respectively, 

/s = c^s/c^i — 1 is the lattice match between the film and substrate along the [100] direction, 
and 03 is the natural lattice constant of the substrate. For p << 1, G as a function of W2, fi, 
fs, and fs can be conveniently used as a state function for the nanorod orientation. When 
G < 0, -Ee«[100] — -EejOOl] < and it is not possible to have nanorods aligned in the [001] 
direction and vice versa. Hence, G = yields a phase boundary that separates the regions in 
the parameter space of {w2, fi, fs, fs) where the vertical or horizontal alignment of nanorods 



is not possible in the c-oriented film on a lattice-mismatched substrate. For a given type of 
impurity doping in the YBCO films, therefore, the tuning of the substrate lattice constant, 
if possible, could result in a transition of the nanorod orientation^^. 

For YBCO epitaxy films on the vicinal substrate, two possible YBCO crystalline config- 
urations near the film/substrate interface have been observed experimentally. In the first 
configuration (see Fig. la), the YBCO ab planes follow the substrate (001) planes on ter- 
races (miscut steps) on the vicinal substrate surface and antiphase boundaries are typically 
formed in YBCO along the c axis at (or near) the step edges of the terraces^"—. In this 
case, the YBCO ab planes can become detwinned with the shorter a axis of the orthorhom- 
bic YBCO layer preferably oriented along the miscut directional^. The percentage of the 
untwinning was observed to increase linearly with the vicinal angle^^ and almost complete 
(70% to 90%) detwinning occurs at 10° vicinal angle in one study^ and less than 1° vicinal 
angle in another study^. In the second observed crystalline configuration (see Fig. lb), 
the YBCO c axis aligns with the normal of the substrate optical surface and the YBCO ab 
planes follow the substrate surface^"—. In this case, very few antiphase boundaries were ob- 
served in the film^ii^^ and a smooth and homogeneous surface morphology was observed by 
the AFM imaging^^. Such a configuration has also been observed in YBCO films on vicinal 
MgO substrates2^. For the doped YBCO films with impurity nanorods on vicinal substrates, 
it is not clear which crystalline configuration dominates near the film/substrate interface, 
although the lack of the antiphase boundary in the samples^^ hints the possibility of the 
second configuration. To study the nanorod configurations, both the crystalline possibilities 
of the YBCO over-layers near the substrate surface were considered in this work. In the first 
possible crystalline configuration, the formation of the antiphase boundaries in the YBCO 
lattice could change the elastic constants as well as the strain field of the YBCO lattice. As 
very few antiphase boundaries have been observed in the YBCO film with impurity nanorod 
inclusions, the effects of the antiphase boundary were not considered here. To model the 
effect of the detwinning of the YBCO ab planes, the lattice constant of the YBCO a axis is 
assumed to decrease linearly from the average lattice constant in the ab plane (a = 3.855 A) 
to its natural value (oqi = 3.823 A) as the vicinal angle increases, i.e. 

0.1 = aoi + (a — aoi)zH{z) (11) 

where 2; = 1 — 0/0^ , 0^ is the vicinal angle at which the ab plane is completely detwinned. 



5 



and H{z) is the Heaviside step function as H{z) = for z < and H{z) = 1 for z > 0. In 
the case that the YBCO c axis is perpendicular to the optical surface of the vicinal substrate, 
the effective lattice constant of the substrate is 

where aos is the natural lattice constant of the substrate and x is a parameter for easily 
switching the direction of the film lattice in the calculation. The film ab planes are parallel 
to the substrate lattice if x = and with the substrate optical surface if x = 1- Including 
both effects of the detwinning and tilting of the film ab planes, the lattice mismatch between 
the film and substrate along the miscut direction ([100] direction) is 

[1 - x(l - cos0)][aoi + (a - aoi)zH{z)] 

In the direction perpendicular to the miscut direction, the mismatch between YBCO and 
STO substrate at the film/substrate interface is negligible. Substituting fs in Eq. (fT3l) into 
Eqs. fl71)- f|T0|) yields G as a function of the vicinal angle </> for the orientation of impurity 
nanorods in the c-oriented YBCO film on the vicinal substrate. Figure 2 plots G = G{(j)) 
for two possible YBCO crystalline configurations on the vicinal substrate. Since G{(f)) > 
when X = 0, the detwinning of the YBCO ab planes is not the root cause of the observed 
transition of the nanorod orientation. In fact, the curves with different 0^ overlaps in Fig. 
2 and, therefore, the change of the strain field in the film due to the detwinning has no 
significant effect on the nanorod orientation. When the YBCO ab planes are parallel to the 
substrate optical surface (x = 1), G{(f)) changes from positive to negative at a threshold 
vicinal angle as shown in Fig. 2. This is when the lattice mismatches between YBCO and 
dopant are substantially altered by the strain in YBCO due to the lattice mismatch with 
the substrate and, consequently, the energetically-preferred nanorod alignment switches from 
the [001] to [100] direction. From Fig. 2, this threshold angle for the onset of the transition 
of the nanorod orientation is found to be (pc — 5.5° for both the BZO and BSO nanorods. 
Experimentally^^, it was observed that the orientation of the nanorods is splayed at = 5° 
and changes into the [100] direction when > 10°, which is in general consistent with 
the model prediction. Since = 5° is very close to 0c, the difference between the elastic 
energies of the horizontal and vertical alignments is very small and the nanorods could 
become splayed. The model calculation has therefore suggested that near the film/substrate 
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interface the film ab planes are parallel with the substrate optical surface and the strain due 
to the lattice mismatch between the film and the tilted substrate lattice is the cause of the 
transition of the nanorod orientation. 

The deformation of the film lattice due to the inclusion of the nanorods can be calculated 
by averaging the principal components of the equilibrium strain over the filn>i^, which can 
be compared with experimental measurement. Because of the different nanorod orientation, 
the deformation of the YBCO lattice is different in the regions of < 0c and > 0c. When 
< 0c, the deformation of YBCO calculated from the equilibrium strain 

ai ~ 2 [A11W2 + (1 + /s + 5/3)Ai2Wi] 2 

(14) 

fcl ^ AiiW2(/3 + Sf-s) Ai3 ^ 

ci ~ 2 [Anit^a + (1 + /a + 5/3)Ai2H 2Ai2 

For BZO or BSO nanorods in YBCO films on STO substrates, + Sf^ > 0. Hence 5ai < 
and 6ci > 0, which represents a compression and expansion of the YBCO lattice along the 
a and c axis, respectively. When > 0c, the deformation of the YBCO lattice is calculated 



Sai 1 



W2{fl + Sfi 



W2 + (1 + /i + 5/3)^1 
Sci I Ai3\ / 6ai \ 



(15) 



Ci \Ai2 J \ ai J 

In this case, the YBCO lattice is expanded along the a axis and compressed along the c 
axis. Note that the ratio of the film lattice deformations along the [001] and [100] direction is 
independent of the properties of the dopant and substrate. Figure 3 plots 6ai/ai and Sci/ci 
as functions of calculated from Eqs. (IHj) and ( IT5l) with fs from Eq. (fT3l) for the case 
of BZO nanorods in the c-oriented YBCO films on vicinal STO substrates. The currently 
available experimental data from the X-ray measurements of the c-axis lattice constant^ was 
also included in the figure. Overall, the model calculation agrees well with the experimental 
data. The discrepancy between the experiment and the model when > 15° is likely due to 
the incoherent growth of YBCO on the vicinal substrate. At a relatively large vicinal angle, 
the lattice mismatch between YBCO and substrate is too big (~ 6.5% at = 20°) to allow 
a coherent film growth. If incoherent lattice growth prevails, the strain in the lattice will be 
released. This is similar to the case of porous YBCO films on vicinal substrates, in which 
a large number of pores appear at approximately 10-15 nm away from the substrate/film 



interface^i^. To further confirm the model prediction on the lattice deformation due to the 
inclusion of impurity nanorods, experimental studies of the lattice deformation in the film 
ab planes are needed. 

In conclusion, a micromechanical model based on the theory of elasticity has been devel- 
oped to study the configurations of the self-organized impurity nanorod array in c-oriented 
YBCO films on vicinal STO substrates. By treating lattice mismatch locally at the inter- 
faces, the strain field due to multiple mismatched lattices of film, dopants, and substrate can 
be considered simultaneously. Including the effect of multiple lattice mismatches is impor- 
tant to the understanding of impurity nanostructures in nanocomposite films when multiple 
lattices are involved. In the case of BZO doped YBCO films on vicinal substrates, two pos- 
sible types of film/substrate interface were studied: perfectly coherent and quasi-coherent. 
The former assumes negligible defects or dislocations on the interface, which results in a sig- 
nificantly strained film lattice. The increase of the lattice strain in the film with the vicinal 
angle of the substrate leads to a transition of the nanorod orientation from the normal to 
the parallel direction of the film. The calculated threshold vicinal angle for the onset of the 
transition for BZO nanorods and the predicated deformation of the YBCO lattice due to the 
inclusion of the nanorods are in very good agreement with the experimental observations. 
This agreement is significant considering no fitting parameter was employed in the model 
and the approach may apply to many other nanocomposites in design and synthesis of novel 
nanostructures. 

This work is supported by NSF and ARO under contract no. NSF-DMR-0803149, NSF- 
DMR-1105986, NSF EPSCoR-0903806, and ARO-W911NF-0910295. 
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Figure 1. Illustration of possible YBCO crystalline configurations near the 
film/substrate interface on a vicinal STO substrate, where the YBCO ah planes par- 
allel to (a) the substrate (001) planes and (b) the substrate optical surface. The 
dashed and vertical solid arrows indicate the [001] direction of the substrate and film, 
respectively. 
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Figure 2. G v.s. vicinal angle for BZO or BSO nanorods in YBCO films on vicinal 
STO substrates. Curves with only difference in 0^ are overlapped. 
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Figure 3. Calculated percentage change of YBCO lattice constants {5L/L) as a func- 
tion of vicinal angle (0) for BZO nanorods in YBCO films on vicinal STO substrates. 
The solid and dashed curve are of the c and a axis, respectively. Open circles are of 
the X-ray measurement of the c-axis. The vertical dashed line marks the threshold 
vicinal angle for the onset of the transition of the nanorod orientation. 
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